Universal Distributions for Growth Processes in 1+1 Dimensions 

and Random Matrices 
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We develop a scaling theory for KPZ growth in one di- 
mension by a detailed study of the polynuclear growth (PNG) 
model. In particular, we identify three universal distributions 
for shape fluctuations and their dependence on the macro- 
scopic shape. These distribution functions are computed us- 
ing the partition function of Gaussian random matrices in a 
cosine potential. 

PACS numbers: 64.60.Ht, 05.40.-a, 68.35. Ct, 81.10.Aj 



Growth processes lead to a rich variety of macroscopic 
patterns and shapes [jjj . As has been recognized for some 
time, growth may also give rise to intriguing statisti- 
cal fluctuations comparable to thermal fluctuations at 
a critical point. One of the most prominent examples is 
the Kardar-Parisi-Zhang (KPZ) universality class In 
essence one models a stable phase which grows into an un- 
stable phase through aggregation, as for example in Eden 
type models where perimeter sites of a given cluster are 
filled up randomly. In real materials, mere aggregation 
is often too simplistic an assumption and one would have 
to take other dynamical modes, such as surface diffusion, 
at the stable/unstable interface into account ||. In our 
letter we remain within the KPZ class. 

From the beginning there has been evidence that in 
one spatial dimension KPZ growth processes are linked 
to exactly soluble models of two-dimensional statistical 
mechanics. Kardar Q mapped growth to the directed 
polymer problem. The replica trick then yields the Bose 
gas with attractive ^-interaction which in one dimen- 
sion can be solved through the Bethe ansatz ||. In ||, 
considerably generalized in Q|, for a particular discrete 
growth model the statistical weights for the local slopes 
were mapped onto the six vertex model. To solve the six 
vertex model one diagonalizes the transfer matrix, again, 
through the Bethe ansatz, which also allows for a study of 
finite size scaling Q . Unfortunately none of these meth- 
ods go beyond what corresponds to the free energy in 
the six vertex model and the associated dynamical scal- 
ing exponent j3 = 1/3. 

In this letter we point out that within the KPZ univer- 
sality class the polynuclear growth (PNG) model plays a 
distinguished role: it maps onto random permutations, 
the height being the length of the longest increasing sub- 
sequence of such a permutation, and thereby onto Gaus- 
sian random matrices p,|lO|. We use these mappings to 



obtain an analytic expression for certain scaling distri- 
butions, which then leads to an understanding of how 
the self-similar height fluctuations depend on the initial 
conditions and to a more refined scaling theory for KPZ 
growth. 

PNG is a simplified model for layer by layer growth 
@. One starts with a perfectly flat crystal in contact 
with its super-saturated vapor. Once in a while a super- 
critical nucleus is formed, which then spreads laterally 
by further attachment of particles at its perimeter sites. 
Such islands coalesce if they are in the same layer and fur- 
ther islands may be nucleated upon already existing ones. 
The PNG model ignores the lateral lattice structure and 
assumes that the islands are circular and spread at con- 
stant speed. The nucleation rate and the lateral speed 
can be set to one by the appropriate choice of space-time 
units. 

We specialize to a one-dimensional surface, returning 
to higher dimensions at the end. The height, h(x,t), at 
time t above the point x on the substrate is counted in 
lattice spacings. The upward steps of h move determinis- 
tically with velocity — f , the downward steps with veloc- 
ity +1, and they annihilate upon touching. Through a 
nucleation event at (x, t), randomly in space-time, h in- 
creases at x by one unit thereby creating a new up-down 
pair of steps. To explain the mapping from PNG to per- 
mutations it is convenient to first use a droplet geome- 
try, where a single island starts spreading from the origin 
and further nucleations take place only above this ground 
layer. The initially flat substrate and other initial condi- 
tions will be handled along the lines of this blueprint. 

We want to compute the height h(x,t) of the droplet. 
Clearly it is determined by the set of nucleation events 
inside the rectangle R( x ,t) = {(#' \t') '■ \ x '\ < t' an d \ x ~ 
x'\ < t - t'}. In lightlike coordinates, r = (*' + x')/V2, 
s = (t' — x')/V2, the rectangle R( x ,t) equals [0, R] x [0, S] 
with R = (t + x)/V%, S = (t - x)/V2. We label the 
nucleation events as (r n ,s n ), n — l,...,N, such that 
0<ri < ■ ■ ■ < riy < R. The corresponding order in the 



second coordinate s, < s 



p(i) 



< 



< s 



p(N) 



< V, defines 



then a permutation p of length N, compare with Fig. [l]. 

There is a simple rule of how to determine the num- 
ber of the layer in which each nucleation event is located. 
Points in layer 1 are obtained by scanning the permuta- 
tion (p(l), ■ ■ ■ ,p(N)) from left to right and marking all 
those entries which are smaller than the so far small- 
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FIG. 1. The height h of a PNG droplet with nucleation 
events corresponding to the permutation (4, 7, 5, 2, 8, 1, 3, 6). 

est. After deleting the subsequence of the first layer, the 
second layer is obtained by repeating this construction. 
One marks those of the remaining entries of the permu- 
tation which are in decreasing order. At the end the 
permutation p has been subdivided into decreasing sub- 
sequences. In the example of Fig. [l]we have the permuta- 
tion (4, 7, 5, 2, 8, 1, 3, 6). The first decreasing subsequence 
is (4, 2, 1), corresponding to the nucleation events in the 
first layer. The remaining subsequences are (7,5,3) and 
(8,6). The height h(x,t) is the number of these subse- 
quences and therefore the length of the longest increasing 
subsequence of p [|l0| . 

In a dual picture one draws a directed path from (0, 0) 
to (x,t), joining nucleation events by straight lines, with 
the restriction that both coordinates r and s are increas- 
ing along the path. Equivalently, the path must be in 
the forward light cone at each nucleation event. This is 
the celebrated directed polymer, cf. for example Q, in 
the context of the PNG model. If to each path we assign 
as negative energy the number of nucleation centers tra- 
versed, then h(x, t) equals the ground state energy of the 
directed polymer. The PNG model is thus in the strong 
coupling regime. 

The nucleation events have density one and are in- 
dependently and uniformly distributed in the rectangle 
R(x.t) with area A = it 2 — x 2 )/2. This induces a Pois- 
son distribution for the length, N, of the permutation, 
ProbjA — n} — e~ x A n /n\, and for a given length n 
each permutation has the same probability, namely 1/n!. 
Thus the problem of computing the distribution of the 
height h(x, t) is converted into determining the statistics 
of the length, I, of a longest increasing subsequence of a 
random permutation. Since to leading order h(x, t) oc t, / 
must be of order vA and the relative fluctuations should 
be of order A 1 / 6 , if we accept (3=1/3 for KPZ growth in 



1 + 1 dimension. 

The same construction can be carried out for an ini- 
tially flat substrate. By translation invariance, it suffices 
to study H{t) = h(0,t). The rectangle R(o,t) is now re- 
placed by the triangle T t = {{x',t') : \x'\ < t-t', t' > 0}. 
To relate to the directed polymer we add the mirror im- 
age relative to t = 0, including the nucleation events, to 
obtain the square R t = {(x',t') : \x'\ < t - \t'\). Then 
2H(t) equals again the ground state energy of the di- 
rected polymer from (0, —t) to (0, t). However the statis- 
tics of nucleation centers inside Rt is constrained to sat- 
isfy the reflection symmetry relative to t = 0. 

For a random permutation with Poisson distributed 
length N, (N) = A, the length I of the longest increasing 
subsequence satisfies the amazing identity 

Prob{/ < to} = e~ x / dUexp (V\Tr(U + U^ 1 ) 

J 771 x m 

(1) 

where the integration is uniformly over all to x to uni- 
tary matrices. A proof can be found for example in 
Jl0|,[ll]]. The partition function in (|l|) appeared before 
in the context of quantum gravity and has a third or- 
der phase transition at to ~ 2\/A with finite size scaling 
governed by the Painleve II equation jl2 13 1. Baik et al 
fTTf prove that I ~ 2\/A + X 1 ^ 6 \2 for large A, where X2 
is a random variable distributed according to the GUE 
Tracy- Widom distribution, i. e. the distribution of the 
largest eigenvalue of a complex hermitian random ma- 
trix @. One has Prob{x 2 < x} = F 2 (x) = e~ B( - x \ 
where g"(x) = u(x) 2 , g(x) — > as x — * oo, and u(x) is 
the global positive solution of the Painleve II equation 
u" = 2u 3 + a;u. Its asymptotics are u(x) ~ W—x/2 for 
x — > — oo and u(x) ~ Ai(x) for x — > oo, Ai(x) the Airy 
function. 

To translate to the PNG model we introduce the 
growth velocity v(u), depending on the macroscopic slope 
u = dh/dx, and the static roughness A(u) Jl5[ , which for 

A(u) = y2 + u 2 in our units. 



PNG are v(u) = V2 + > 
Then 

h{v'{u)t, t) ~ (v(u) - uv'{u))t + {\v"{u)A(u) 2 i) 1/3 X2 

(2) 

in the limit of large t. We emphasize that all nonuniversal 
factors are given through the model dependent quantities 
v(u), A{u) and remark that ( B) i s also confirmed by the 
rigorous result of Johansson |16[ | for a discrete growth 
model equivalent to the totally asymmetric simple exclu- 
sion process. 

For the flat substrate one might expect to have the 
same fluctuation law as for the droplet, since in both 
cases the mean curvature vanishes on a microscopic scale. 
A result of Baik and Rains |l7j tells us however that the 
fluctuations are GOE [14| . More precisely, there is a sim- 
ilar formula for Prob{Z < to} as (|lj) in the case that the 
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random permutation p is reflection symmetric relative 
to the anti-diagonal, p(N + 1 - p(k)) = N + 1 - k. The 
asymptotic analysis of |l7j results in I ~ 4\/A+(2A) 1 / 6 Xi, 
where xi is distributed as the largest eigenvalue of a real 
symmetric random matrix. Translated to the surface this 
means 

H{t) = V2t + (t/V2)V 3 Xl . (3) 

One has Prob{2 2 / 3 xi < x} = F 1 (2- 2 / 3 x) = 
e -(/(a ! )+fl(x))/2 ) g ( x j as above and ft( x j = _ u ( a; ) ) /(a;) _> 

for a; — ■> oo. The distributions of \2 and X i are plotted 
in Fig. ||. Superimposed are Monte Carlo data for the 
PNG model, which differ distinguishably from the ana- 
lytical curves only at the tails where statistics becomes 
bad. We conclude that the droplet and the flat substrate 
have the same scaling form but distinct universal distri- 
butions. 

The flat substrate, although used in many simulations, 
and the droplet are rather special as initial conditions. 
From a statistical mechanics point of view stationary 
growth would be regarded as singled out, which for PNG 
corresponds to initial conditions where the up and down 
steps are random with densities \[2 each. Physically, 
another natural initial condition is to have a staircase 
configuration representing a tilted surface. In addition 
we could have sources, for example additional nucleation 
events at the origin. The mapping to the directed poly- 
mer works as before. Our crucial observation is that such 
other initial conditions translate in essence to defect lines 
and/or boundary potentials for the directed polymer. 

To illustrate, we discuss only one special geometry. As 
for the droplet we consider random nucleations of density 
one in the square R(o,t)- I n addition there are random 
nucleations at the two lower edges {s = 0} and {r = 0} 
with constant line densities p+, resp. /?_. Thus the path 
of minimal energy, with starting point at (0,0), sticks 
for a while at one of the two edges and then enters the 
bulk to reach (0, t) eventually. If p+ < 1, p- < 1, it 
does not pay to stay at the edges, and from the bulk we 
have GUE energy fluctuations according to (||). On the 
other hand if p = max{p + ,p_} > I, the optimal path 
stays for a length t(l — 1/ ' p 2 )j\[2 at the edge with the 
higher density. Since the edge events are random, the 
i 1 / 3 bulk fluctuations are dominated by the Gaussian \fi 
edge fluctuations. Parenthetically we remark that for 
regularly spaced edge points one should recover GOE. 

The length distribution along the critical lines p+ — I , 
P- < 1, resp. p+ < 1, p- = 1, was identified by Baik 
and Rains |l8) , in a generalization of the techniques in 
[ pTf . They obtain GOE 2 fluctuations, i.e. the distribu- 
tion of the maximum of two independent GOE random 
variables. The path of minimal energy stays for a length 
of order t 1 / 3 at the density one edge. At the critical point 
p+ = p- = 1, the polymer has a choice between the 
left and right edge. By a limiting procedure one obtains 
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FIG. 2. From left to right: the probability densities of the 
universal distributions X2, Xii an d Xo f° r curved, flat, and 
stationary self-similar growth, respectively. 

pl| the universal distribution for the energy fluctuations, 
Fo(x) — Probjxo < x}, with 

F (x) = [l-(x + 2f" + 2g")g']e-^\ (4) 

An interpretation in terms of the eigenvalue distribution 
of random matrices has yet to be found. In Fig. we plot 
the distribution of xo- Superimposed are simulation data 
for the PNG model, taken before the analytic result had 
been obtained. The first four moments of Xh i = 0> 1> ^' 
are listed in Table |[ Of interest are also the asymp- 
totics of the probability densities Fj(x). From Painleve 
II we obtain — \ogFj(x) — Cj\x\ 3 /12 for x — > — oo and 
— \ogFj(x) — djX 3 / 2 /3 for x — » oo up to logarithmic cor- 
rections with prefactors Cj = 1,2,1 and dj = 2,4,4 for 
j = 0, 1, 2, respectively. 

We have to translate back to surface growth. For sta- 
tionary growth with zero slope, in the space-time picture, 
the height lines cross the forward light cone with the den- 
sities p + = 1 = p_ and the intersection points are Poisson 
distributed |2(| . Thus for the directed polymer with edge 
densities the critical point is precisely stationary growth 
with zero slope. If p+p~ = I , p+ ^ I, we have also sta- 
tionary growth but now with slope u = (p_ — p + )/\/2. 
As argued already the fluctuations along the line x = 
are then Gaussian ^ft. For the t 1 ^ 3 fluctuations one has 
to record height differences along the line {x = v'(u)t\, 
as can be seen from a similarity transformation. In Fig. y 



TABLE I. Mean, variance, skewness, and kurtosis for the 
distributions of X2, Xi> an d Xo as determined by numerically 
solving Painleve II [19]. {x n )c denotes the n'th cumulant. 





curved (X2) 


flat ( X i) 


stationary (xo) 


<x> 


-1.77109 


-0.76007 





0.81320 


0.63805 


1.15039 


(x%Kx 2 )T 


0.2241 


0.2935 


0.35941 


0.09345 


0.1652 


0.28916 
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FIG. 3. Droplet with boundary sources of intensity p_ < 1 
and p+ = 1 (critical). The dashed line corresponds to the free 
droplet. The shaded region is the extra mass due to sources. 

we illustrate the macroscopic shape for general boundary 
sources. If = = p_, we have the droplet discussed 
before. Nonzero boundary sources enforce flat segments 
tangential to the droplet shape. The profile at x = is 
curved for p + < 1, p_ < 1, flat otherwise, the marginal 
case corresponding to the critical lines. 

Our detailed study of the PNG model suggests the fol- 
lowing scaling theory for all growth models in the KPZ 
universality class. First of all we require a self-similar 
macroscopic shape. Locally this leaves only two possibil- 
ities, either a flat piece or a curved piece with a shape 
determined through the slope dependent growth velocity 
pl| . We draw a ray from the center of symmetry. If 
the surface at the point of intersection with the ray has 
non-zero curvature, then the height fluctuations in this 
direction are GUE with scaling form (||). If the curva- 
ture is zero, we have to know the roughness of the ini- 
tial conditions, i. e. \h(x, 0) — h(0, 0)| oc \x\ a with rough- 
ness exponent a, and/or the corresponding roughness for 
boundary sources. If a = the height fluctuations are 
GOE and the general scaling form is as in (Q) with \i 
replaced by xi- If a = 1/2 the height fluctuations are 
Gaussian with variance proportional to t, except along 
the line {x = v'(u)i\, where they again have the scaling 
form <M) with the random variable %2 replaced by Xo, as 



defined in 



The intermediate cases < a < i have 



not been studied systematically. Also the fluctuations 
at the endpoints of flat pieces have still to be classified. 
There are two exceptions. One is the case of Fig. ||, 
which has GOE 2 and the second one is the half-droplet 
with an external source at x = 0. Translating JItJ to 
PNG one finds GOE, GSE, and Gaussian depending on 
the strength of the source. 

Our constructions carry over immediately to higher di- 
mensions, as can be seen most directly in the polymer 
picture. The square is replaced by a (d+ l)-dimensional 
(hyper)cube with uniformly distributed nucleation cen- 
ters, the polymer running from the lower to the upper tip. 
For the PNG model this corresponds to droplet growth 
with islands having the shape of a regular simplex (a tri- 



angle in 2 + 1, a tetrahedron in 3 + 1, and so on). One 
axis of the cube defines the order 1, . . . , N, while the re- 
maining d axes define the permutations Pi(l), ■ . . ,Pi(N), 
i = 1 , . . . , d. Increasing means now increasing in all co- 
ordinates, i.e. j < j' and Pi(j) < Pi{j') for all 1 < i < d. 
The length of the longest increasing subsequence equals, 
again, the height of the droplet. At present we study 
numerically the statistics of this length with the goal to 
have information on scaling more precise than the one of 
previous investigations Q . 

In conclusion, we have obtained distinct scaling func- 
tions for the PNG model, which depend on the choice of 
initial conditions. By universality we argued that from 
the knowledge of the self-similar curvature one can infer 
the type of height fluctuations. It would be of interest 
to study also joint probability distributions of the height 
at distinct space-time points. Perhaps such a program 
could identify the universal field theory hiding behind 
KPZ growth in one dimension. 

We are grateful to J. Baik and E. M. Rains for making 
their results Jl^] available to us prior to publication and 
to C. Tracy for help on Painleve II. 
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